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1. Introduction 

Gravitational lenses are becoming a precious mean of probing the matter 
distribution in the Universe. From the search of matter in the Galactic halo 
to the study of the large-scale structures of the Universe, the gravitational 
lens effects offer a unique alternative to light surveys and are now widely 
used. This evolution is due in particular to the use of new observation 
devices, such as the wide field CCD cameras. 

The aim of this course is to study the effects of gravitational lenses in 
those different astrophysical contexts. These notes are voluntarily focused 
on the fundamental mechanisms and the basic concepts that are useful to 
describe these effects. The observational consequences will be presented in 
more details in these proceedings by Y. Mellier (see also his review paper 
1998). Related textbooks are, 

— "Gravitation" by Misner, Thorne and Wheeler (1973) for General Rel- 
ativity and in particular for the presentation of the geometric optics. 

— "Large-scale Structures of the Universe" by Peebles (1993) for the de- 
scription of the large-scale structures of the Universe. 

— "Gravitational lenses" by Schneider, Ehlers &: Falco (1992) for a gen- 
eral (but rather mathematical) exhaustive presentation of the lens 
physics. 

The content of these notes is the following. In the first section I describe 
of the basic mechanisms of gravitational lenses, techniques and approxima- 
tions that are usually employed. The second section is devoted to the case 
of a very simple deflector, a point-like mass distribution. This corresponds 
to microlensing events in which the deflectors are compact objects of a frac- 
tion of a solar mass that may populate the halo of our Galaxy. The last two 
sections are devoted to cosmological applications. After a presentation of 
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Figure 1. The geometrical relationship between the deflection angle 6 and the displace- 
ment angle /J. 

the geometrical quantities that are specific to cosmology, I will present the 
various phenomena that can be observed in this context. Finally I describe 
the weak lensing regime. This is a rapidly developing area that should even- 
tually allow us to map the mass distribution in the Universe. I explore how 
this can be used to constrain the cosmological parameters. 

2. Physical mechanisms 

The physical mechanisms of gravitational lenses are well known since the 
foundation of General Relativity. Any mass concentration is going to deflect 
photons that are passing by with a fraction angle per unit length, S6/Ss, 
given by 



where the spatial derivative is taken in a plane that is orthogonal to the 
photon trajectory and cf) is the Newtonian potential^. 

^We will see in section 4 what is its meaning in a cosmological context. 
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2.1. BORN APPROXIMATION AND THIN LENS APPROXIMATION 

In practice, the total deflection angle is at most about an arcmin. This is the 
case for the most massive galaxy clusters. It implies that in the subsequent 
calculations it is possible to ignore the bending of the trajectories and 
calculate the lens effects as if the trajectories were straight lines. This is 
the Born approximation. 

Eventually, one can do another approximation by noting that in gen- 
eral the deflection takes place along a very small fraction of the trajectory 
between the sources and the observer. One can then assume that the lens 
effect is instantaneous and is produced through the crossing of a plane, the 
lens plane. This is the thin lens approximation. 



2.2. THE INDUCED DISPLACEMENT 

The direct consequence of this bending is a displacement of the apparent 
position of the background objects. This apparent displacement depends 
on the distance of the source plane, Dqs^ and on the distance between the 
lens plane and the source plane Dls- More precisely we have (see Fig. 1), 

/3 = a--^— — ^ — Va ( [ ds (p{s,a)) (2) 



c2 DosDoL 

where a. is the position in the image plane, (3^ is the position in the source 
plane. The gradient is taken here with respect to the angular position (this 
is why a Dql factor appears). The total deflection is obtained by an inte- 
gration along the line of sight, assuming the lens is thin. In a cosmological 
context the exact expressions of the angular distances are not trivial, they 
depend on the local curvature of the background. 

3. The case of a point-like mass distribution 

3.L MULTIPLE IMAGES AND DISPLACEMENT FIELD 
The potential of a point-like mass distribution is given by, 

for an object of mass M. Let me calculate the instantaneous deflection 

angle at an apparent distance r. We suppose that the impact parameter of 
the trajectory is r and x is the abscissa to the point of the trajectory that 
is the closest to the lens (see Fig. 1). Along the trajectory the potential is 
given by, 

^/ A -GM ... 
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Figure 2. Position of the Einstein ring for a point-like mass distribution. 



Then the deflecting angle is given by, 

- o GMr 
6x c2(r2 + a;2)3/2- ^ ^ 

The total deflection angle 9 is given by the result of the integration of this 
quantity with respect to x. It gives, 

6 = (6) 

This is a well known expression which served as the first test of General 
Relativity with the observed displacement of stars around the solar disc 
(Dyson et al. 1919). 

It implies that the true position of an object on the sky, /3, is related to 
its apparent position, a, with 



p = d-^d (7) 
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where (5 and a are 2D angular position vectors and Re the Einstein 
radius, 



One can see that when = R\ the lens and the background objects are 
necessarily aligned. It implies that, since the optical bench is symmetric 
around its axis, the observed object appears as a perfect ring (see Fig. 2). 
It is worth noting that for this potential, except for this particular position, 
all background objects have two images. This is however quite specific to a 
point like mass distribution which has a singular gravitational potential. 

The problem is that none of these features are observable when the lens 
is a star. Let for example assume that we have a one solar mass star in the 
halo of our galaxy (therefore at a distance of about 30 kpc). The apparent 
size of such a star is about 10~^ arcsec. Its Einstein ring is about 10~^ 
arcsec^. None of these dimensions are accessible to the observations (the 
angular resolution of telescope is at best a few tens of arcsec). The Einstein 
radius is therefore much too small to be actually seen! 

Note however that these numbers show that the point-like approxima- 
tion is entirely justified for a star (the Einstein ring is much more bigger 
that the apparent size of a star) . Simple examination of the scaling in those 
relations shows that this would not be true for massive astrophysical objects 
such as galaxies or galaxy clusters. 

The detection of gravitational effects due to stars should then be done 
by another mean: the amplification effect. 

3.2. THE AMPLIFICATION MATRIX 

The case of circular lenses has already given us a clue: when the source is 
precisely aligned with the lens, the image is no more a point but a circle. 
One consequence is that the observed total luminosity is much larger than 
what would have been observed without lenses. The effect is basically due 
to the variations of the displacement field with respect to the apparent po- 
sition. These variations induce a change of both the size and shape of the 
background objects. To quantify this effect one can compute the amplifica- 
tion matrix A which describes the linear change between the source plane 
and the image plane. 



^To do this calculation it is useful to know that the horizon of a one solar mass black 
hole, r = 2GMq/c^, is about 3 km. 
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Its inverse, A~^, is actually directly calculable in terms of the gravitational 
potential. It is given by the derivatives of the displacement with respect to 
the apparent position, 



- da, ~ 



Dls 



*j Dos DoL 

In case of a point-like mass distribution it is easy to see that, 



1 - 



a? 



(10) 



(11) 



The amplification effect for each image is given by the inverse of the deter- 
minant of the amplification matrix computed at the apparent position of 
the image. The amplification factor is usually noted /x, 



/x = l/det(^-^). 



(12) 



In case of the point-like distribution we have (the calculation is simple at 
the position ai = a, 0:2 = 0), 



a 



(13) 



for each image. The total amplification effect is given by the summation of 
the two effects for the 2 images, 



(14) 



where h is the impact parameter of the background object in the source 
plane. The amplification effect is obviously dependent on the impact pa- 
rameter. If it is changing with time, this effect is detectable. 



3.3. THE MICROLENSING EXPERIMENTS 

The microlensing experiments are based on this effect. When a compact 
object of the halo of our galaxy reaches, because of its proper motion, the 
vicinity of the light path of a background star (from the SMC or the LMC) 
the impact parameter is changing with time and can be small enough to 
induce a detectable amplification (when u is about unity, the amplification 
is about 30%). In practice one observes changes in the magnitude of the 
remote stars that obey specific properties, 

— the time dependence of the amplification is symmetric and has a spe- 
cific shape; 
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— the amplification effect is unique; 

— the magnitude of the amphfication effect is the same in all wavelengths. 

The time scale of such an event is about a few days to a few month 
depending on the mass of the deflectors. Currently a fair number of such 
events have been recorded (see contribution of J. Rich, these proceedings) 
and constraints on the content of our halo with low massive compact objects 
have been put. 



4. Gravitational lenses in Cosmology 



The extension of the lens equations to a cosmological context raises some 
technical difficulties because the background in which the objects are em- 
bedded is not fiat. The aim of this section is to clarify these points. However, 
readers that are not familiar with cosmology can jump to section 5. 

The basic equations, that describe jointly the evolution of the expansion 
parameter and the mean density, are the following. 



3 d 



-AttG pa + Aa; 



8ttG 



pa^ — ka + 



Aa^ 



(15) 
(16) 



where A is a possible cosmological constant and k a possible curvature term. 
The Hubble constant reads, 



(17) 




To simplify the discussions the reduced quantities are introduced, 

s i = :r-r:^ p and A = 



3^2^ — 3^2- 
They have an index when they are taken at present time. 



(18) 



4.1. THE ANGULAR DISTANCES 

We consider an object of size I (either because of its proper size or because 
of a peculiar physical process) at redshift z. When this size I is seen under 
an angle a, then by definition, 

« = (19) 

where Pq is the angular distance. This is the distance at which this object 
would be in an Euclidean metric. What is then the relationship between 
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I?o and z? The Friedmann-Robertson- Walker metric is given by, 

ds^ = -c^ di^ + a\t) (-^^ + x^^e" + sin^ 0d<^2^ . (20) 

\ -L K % j 

The fact that the size of this object is / means that it takes a time interval 
dt = l/c for light to travel to one end to the other. The corresponding angle 
a = AO can be obtained by writing ds = with Ax = 0, d</? = 0, which gives 

l = a{t)xa. (21) 

The angular comoving distance is thus given by ^^^{z) = x{z). The expres- 
sion of x{z) can then be computed by the relation ds = along the line of 
sight with d^ = and dip = 0, 

/■*! cdi _ dx ^ ^^^^ ^^^^ 

Jto a Jo Vl — kx'^ 

For an open Universe, A; < 0, and we have, 

x{z) = Vo{z) = sinh[V^ x(^)]- (23) 

Obviously when A; = 0, x = x- 

The relation x(z) depends on the function a{t) for a given cosmology, 
and therefore on the matter content of the Universe, on A and k . For 
instance, for an Einstein-de Sitter Universe, in the matter dominated era, 
we have 

2/3 



>m = {^) , (24) 



so that, 

x{z) = r t-y^dt = 3{tl/' - tl/')tf', Ho = ^, (25) 
Jto ^ to 

and eventually. 

This is the comoving angular distance for an Einstein-de Sitter Universe. 
More generally we have an explicit solution if A = only. Finally the lens 
equation also requires the angular distance between two different redshifts 
zi and Z2- The calculation is actually quite simple. This distance is given 
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formally by 'Dq{z2) when it is calculated at a time when the observer is at 
redshift zi. k being time independent we should formally have, 



1 



'Do{zi,Z2) = ^=sinh[V^ (x(^2,^i))], 



/-k 



and to compute x(-22, zi) one only needs to remark, 

r*2 cdt 



X{z2, zi)= f — = x{z2) - x{zi), 

Jti a 



(27) 



(28) 



which gives the expression of the angular distance we need. Eventually, it 
is fruitful to notice that, 

cdt cdz 



H ' 



which gives, 

'^0{Z1,Z2) 

and 

-Doiz) = 



sinh 



J Z\ 



^2 dz 



H{z) 



sinh 



Hq ^1 — i^o — Xq 



dz' 




(29) 



(30) 



(31) 



The whole geometrical part of the lens equation is thus established. 



4.2. GEOMETRIC OPTICS IN A WEAKLY INHOMOGENEOUS UNIVERSE 

What is now the source term for the deflection angle? We should first no- 
tice that in absence of lenses the light rays follow the geodesies of the 
Priedmann-Robertson- Walker metric. And in the applications we are inter- 
ested in, the metric fluctuations arc always weak. These fluctuations are 
given by, G M/{Rc^). For instance, for 

- 1 star: M = IMq, R = 710^km, 6(j) « 10"^; 

- 1 galaxy cluster: M = 10^^ Mq, R= IMpc = 3 lO^'^km, 6(f) ^ 10'^. 

The metric inhomogeneities are thus always extremely weak, even in the 
most extreme cosmological situations. 

Following^ Sachs (1961), we consider two nearby geodesies, C and C, 
in a light bundle in an FRW Universe with small metric fluctuations. We 
denote the bi-dimensional angular distance between C and C' as it is 
seen by the observer. This is the distance in the image plane, that is the 

^See Misner, Thome and Wheeler for an exhaustive presentation of the geometric 
optics. 
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Figure 3. Sketch showing the geometrical quantities that appear in Eq. 32 

difference between tlie angular coordinates with whicli the photons arrive. 
We denote ^i{z) the real distance between C and C at redshift z (see Fig. 
3). It implies that the geodesies are straight enough so that light always 
travels towards the observer. We also assume that the deflections are small 
enough so that it is possible to make the small angle approximation, 



that is that we assume the the position vector can be obtained by a 
simple linear transform of the angular coordinates. For an homogeneous 
space Vij{z) is simply given by Vq{z) where 6^ is the Kronecker symbol. 
Obviously Vij changes as a function of redshift along the trajectories. The 
"virtual" angular position in the source plane is then given by the ratio 
of the real distance (at time of light emission for instance) by the angular 
distance of the emitter in an homogeneous space, 



The amplification matrix, or rather its inverse, A , is then given by. 




(32) 




(33) 




(34) 
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for a source plane at redshift z. 

Sachs (1961) gave the master equation which governs the evolution of 
the distance between the geodesies. The derivation of this equation goes 
beyond these lecture notes and I give only the final answer (as it has been 
given by Seitz, Schneider and Ehlers 1994), 

d^[ay z)] ^ ^) .(^^ ^) ^35) 

where the derivatives are taken with respect to r/, 

d?7 = --^ = -adi, 7?(z = 0) = 0. (36) 
H(a) 

with the boundary conditions, 

The matrix T^jj represents the tidal effects. It can by written in terms of 
the gravitational potential given by, 

A^^ = 47r G p Jmass • (38) 

The Laplacian is taken with respect of the comoving angular distances. We 
have 

Hla^\Q 1 J H^a^ [ <j>,2i 0,22 ) " ^ ^ 

Since SirGpa^ = 3HqQq, (e.g. Eq. 18) for an homogeneous Universe we 
have, 

n^^ = -l{i + zfno6i^ (40) 

(the superscript (0) means here that it is the value of TZ without perturba- 
tions). In this case the matrix Vij is proportional to S-j and we have, 

^!tLpM._|(i + ,)4„„„„„. (41) 

We recover in fact the comoving angular distance the expression of which 
we know, 

c r dz' 1 
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with 

E{z) = ^ = JXo + (1 + ^)2(1 -no- Ao) + (1 + z)3no. (43) 

This integral has a closed form only when the cosmological constant, A, is 
zero. 

4.3. THE LINEARIZED EQUATION OF GEOMETRIC OPTICS 

We can remark that the equation (35) is not linear since V is not simply 
proportional to TZ. This expresses the fact that the deformation of the an- 
gular distance is made all along the light trajectory by multiple deflections. 
The general resolution of Eq. (35) is in general very complicated. It can 
however handled when it is linearized. Let's assume we can expand Vij 
with respect of the local density contrast, 

Vij{z)=Vo + vff + ... (44) 

It implies that, at first order, 

d'^[aVi]\p,z)] .n^ ^ ^ 

- aiz)n^\p,z)vi'JiP,z) = -30o(l + zrVo{z)^,j{P,z) 

(45) 

with. 



z=0 



and we define the field (p so that, 

^xV^ = ^mass(/3, Z) = . n - (47) 

47r G p 

To solve this differential equation it is easier to write it with the variable 
z. It then reads, 

d^[aVif0,z)] ^ 1 dEiz)dvif0,z) 
dz^ E{z) dz dz 

l + zE{z) dz '^*^^^'^^ + 2 E^{z) "^'^ ^^'^^ 

-5 ^f,^^^'^ Vo0z)<p,ij0,z). (48) 
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The differential liomogeneous equation it is associated with has two known 
solutions. One describes the angular distance Vq, the other is given by, 



Uoiz) = 



cosh 



Vi-no- Ao /' 
Jo 



dz' 



rio — Ao 

The general solution of 48 then reads, 



E{z') 



(49) 



V\f0,z) = -3noj^'dz' 



Uoiz')Voiz) 



EHz) 



U(,{z')Vo{z')-Uoiz')V'^{z')' 

(50) 

which, after elementary mathematical transforms, gives. 



dz' 



X sinh 



dz" ' 



Jz 



E{z") 



{l + z')VQ{z')ip,ij{z'). 



(51) 



It can be rewritten by introducing the physical distance % along the line of 
sight. We eventually have. 



A-^{z) = ld- 



3Qo fx{z) ^Vo{z',z)Vo{z') 



dx' 



-Doiz) 



{l + z')<p,,j{z'), (52) 



where the angular distances 'Do{z) and "Dq^z, z') arc comoving. This equa- 
tion actually gives the expression of the amplification matrix for a non- 
trivial background. We find that the amplification matrix is given by the 
superposition of lens effects of the different mass layers. Wc can remark 
that the lens term is given by the gravitational potential, 0, that is by the 
potential the source term of which is given by the density contrast. 

Note finally that this equation is valid in two limit cases, either for a 
single lens plane with an arbitrary strength or the superposition of any 
number of weak lenses. This equation naturally extents the previous result, 
(10), obtained for a single lens in an Euclidean background. The higher 
orders of Eq. (35) give the intrinsic lens coupling effects (i.e. their nonlinear 
parts). We will not consider them here. 



5. Galaxy clusters as gravitational lenses 

The study of galaxy clusters has become a very active field since the dis- 
covery of the first gravitational arc by Soucail et al. (1988) in Abell cluster 
A370. Galaxy clusters give the most dramatic example of gravitational lens 
effects in a cosmological context. The difficulty is however to describe the 
shape of their mass distribution. 
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5.1. THE ISOTHERMAL PROFILE 

For an isothermal profile we assume that the local density p{r) behaves 
like, 

pW=Po(^)"'. (53) 

With such a density profile the total mass is not finite. So this is not a 
realistic description but it is a good starting point for the central part of 
clusters. It is actually more convenient to parameterize the depth of a po- 
tential well with the velocity dispersion it induces. The velocity dispersion 
is due to the random velocity that particles acquire when they reach a sort 
of thermal equilibrium. Such a dispersion is in principle measurable with 
the observed galaxy velocities along the line of sight. The velocity disper- 
sion is related to the mass M(< r) of the potential well that is included 
within a radius r, 

2, , GM{< r) 
a^[r) ^ ^ -. (54) 

In case of a isothermal profile, the velocity dispersion is independent of the 
radius and we have 

a'^ = 2tiG porl. (55) 
The integrated potential along the line of sight is given by, 

ip(r) = 2ir r. (56) 

As a consequence the amplitude of the displacement is independent of the 
distance to the cluster center and 



- ^ 47r Dls 2 " /(--x 
13 = a- — — — cr -. (57) 
Dos ct 



The position of the Einstein ring is obviously given by, 

47r Dls 2 frQ\ 
a = Re = — — o- , (58) 
Dos 

which depends both on the velocity dispersion and on the angular distances. 
The number of images depends in this case on the value of the impact pa- 
rameter. If it is too large (i.e. larger than Re) then each background object 
has only one image. For a galaxy cluster of a typical velocity dispersion of 
500 km/s, and for a source plane situated at twice the distance of the lens, 
the size of the Einstein ring is about 0.5 arcmin. It is interesting to note 
that the size of the Einstein ring is directly proportional to the square of 
velocity dispersion (in units of c^) and to the ratio Dls/Dqs- 
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Figure 4- Graphical determination of the position and number of images from the shape 
of the potential. 



The amplification matrix reads, 



where we have, 



As a result the amplification is given by, 

H = -^. (61) 

1 — X 

Once again the amplification becomes infinite when x — > 1 that is, close to 
the critical line. 



5.2. THE CRITICAL LINES FOR A SPHERICALLY SYMMETRIC MASS 
DISTRIBUTION 

The two previous cases correspond to specific profiles. In this part 1 only 
assume a spherical symmetric profile. The displacement is then given by 
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the derivative of the potential, 

^ _ dip d 

P = d--f-. 62 

dr a 

It is interesting to visuahze this relation with a graphic representation. This 
is proposed in Fig. 4. The number and position of the images of a given 
background object are given by the number of intersection points between 
the curve and a straight line of slope unity. This is a direct consequence of 
the relation, 

b-a = -^ (63) 
da 

when the potential is computed along a given axis that crosses the cluster 
through the center and a and b are the abscissa on this axis of one given 
object in respectively the source and the image plane. 

The interesting quantity is also the amplification matrix that indicates 
the position of the critical lines. In general this matrix reads. 



1 - ^ 

or-' 

1-^^ 



when it is written in the basis (ev, eg). Then the amplification is infinite in 
two cases, when 

^ = 1 or -— = 1 (65) 

dr"^ r dr 

The second eigenvalue corresponds to the same case as for a singular isother- 
mal profile. At this particular position the source forms an Einstein ring. 
The first eigenvalue, however, is associated with an eigenvector that is along 
the X direction, that is along the radial direction. It means that the "arc" 
which is thus formed is radial. It graphically corresponds to the case of two 
merging roots. It is therefore directly associated with the behavior of the 
potential near the origin. 

5.3. THE ISOTHERMAL PROFILE WITH A CORE RADIUS 

Let us consider a simple case where the projected potential is made regular 
near the origin, 

ifir) = ifo ^1 + (r/re)2 (66) 
The constant ipo is related to the velocity dispersion with 



(^0 = — p; (67) 

Dos Tc 
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Figure 5. Shape of the caustic Unes (thick hrics) and critical hnes (thin lines) for an 
elliptic potential and for different values of the central potential. 

where a is here the velocity dispersion at a radius much larger that Tc (the 
velocity dispersion decreases to zero at the origin in this model). This is a 
more realistic case. It is interesting to note that in this case the potential is 
not necessarily critical (there may be no region of multiple images, see Fig. 
4) . When it is critical the discovery of a radial arc is an extremely precious 
indication for the value of the core radius. 



5.4. CRITICAL LINES AND CAUSTICS IN REALISTIC MASS 
DISTRIBUTIONS 

In realistic reconstructions of lens potential however, it is very rare that 

the lens is circular. Most of the time the mass distribution of the lens is 
much more complicated. It induces complex features and series of multiple 
images. 

The simplest assumption beyond the spherically symmetric models is 
to introduce an ellipticity e in the mass distribution (Kassiola &: Kovner 
1993), 

(p = (po^l + rlJrl with rl^ = 7^332 + ^ +e)2 " ^^^^ 
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To understand the physics it induces one should introduce the caustics 
and critical lines. The critical lines are the location on the image plane of 
the points of infinite magnification. The caustics arc the location of these 
points on the source plane. These points are determined by the lines on 
which det(^^^) = 0. It means that arcs are along the critical lines and 
that they are prodiiccd by background galaxies that happen to be located 
on the caustics. On Fig. 5 one can see the shapes of the critical lines and 
caustics for different depth of the potential (68) (or equivalently for different 
positions of the source plane). 

Eventually the reconstruction of galaxy cluster mass maps requires the 
use of more complicated models and it can be necessary to perform non- 
parametric mass reconstructions. Recent results have been obtained by Ab- 
delSalam et al. (1997) for few clusters. 

6. The weak lensing regime 

In this section I consider the possibility of using the lens effects to probe 
the large-scale structures of the Universe. The difficulty is here that the 
distortion induced by the lenses can be very small. The projected potential 
should then be reconstructed with a statistical analysis on the deformation 
measured on a lot of background objects. Let me define more precisely the 
different regimes. 

6.1. THE MATHEMATICAL DESCRIPTION OF THE WEAK LENSING 
REGIME 

Depending on the magnitude of the lens effect, different regimes are possi- 
ble: 

— The strong regime is such that several light path are possible between 
the sources and the observer. It thus induces multiple images, and the 
images of background galaxies are often extremely distorted: this is 
the regime of the giant arcs. The one that has been investigated so far. 

— The regime of arclets correspond to a case where there are no multiple 
images although a significant distortion of the background objects can 
be observed. 

— The weak lensing regime corresponds to cases of modest distortion 
(typically a few %). Such effect cannot be detected with a single object 
and therefore should be measured in a statistical way, by averages over 
a large number of background galaxies. 

In all cases the displacement field is not directly observable. In the weak 
lensing regime, the deformation only is measurable. For slightly extended 
objects such as background galaxies the deformation in shape is induced by 
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the variations of the displacement field. It can actually be described by the 
amplification matrix, the components of its inverse are in general written, 

taking advantage of the fact that it is a symmetric matrix. The components 
of this matrix are expressed in terms of the convergence, k, (a scalar field) 
and the shear, 7 (a pseudo vector field) with 

K=^V'^tP; 71 = ^(V',!! - V',22) ; 72 = V',12> with V = 2— -^^^^ — (f). 
2 2 ^JOS Uql 

(70) 

The convergence describes the linear change of size and the shear describes 
the deformation. The consequences of such a transform can be decomposed 

in two aspects: 

— The magnification effect. Lenses induce a change of size of the objects. 
As the surface brightness is not changed by this effect, the change of 
surface induces a direct magnification effect, ^. This magnification is 
directly related to the determinant of A so that, 

IX = det{A) = 1/ [(1 - - 7^] . (71) 

— The distortion effect. Lenses also induce a change of shape of the back- 
ground objects. The eigenvalues of the matrix determine the di- 
rection and amplitude of such a deformation. 



6.2. THE MAGNIFICATION EFFECT 

In the weak lensing regime, the observational consequences of the magnifi- 
cation effect is a combination of a change the apparent area of the objects, 
that makes their detection easier, and their dilution (the total area is en- 
larged as well). The mean local number density of galaxies is then related 
to the slope of the galaxy counts, a, through 

n(7) =n„//2-5"-i with a = ^^?^, (72) 

dm 

where m is the apparent magnitude (in a given band) and Ug is the mean 
number density of galaxies in the absence of lenses. Whether the number 
density of galaxies increases or decreases in magnified area thus depends 
on a and consequently on the selected population of objects. 

Such an effect has been advocated (Broadhurst 1995 and Broadhurst 
et al. 1995) as a way to detect the lens effect. In general this method 
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suffers from the fact that the background galaxies have intrinsic number 
density fluctuations. It is therefore more appropriate for mapping the mass 
fluctuations in galaxy clusters where the magnification effect is large enough 
to dominate. In galaxy clusters it is a cheap way to map the mass profiles. 
More sophisticated analysis allow even to have access to the cosmological 
constant A by probing the extension of the depletion area (Fort et al. 1997, 
see contribution of Y. Mellier) . 

6.3. THE GALAXY SHAPE MATRICES TO MEASURE THE DISTORTION 
FIELD 

The distortion effects change the shape of the background objects. The 
objects appear elongated along the eigenvalues of the amplification matrix. 
When background objects are only moderately extended (this excludes the 
case of arcs), their shapes can be described by the matrix, 

s = J d^eeiOjiie). (73) 

It is easy to relate the shape matrix in the source plane to the one in the 
image plane. This is obtained by a simple change of variable that uses the 
fact that the surface brightness of the objects is not changed. It implies 

5^ = ^-1.5^.^-1. (74) 

By averaging over the shape matrices in the source plane, assuming the 
intrinsic shape fluctuati ons are not correlated, one can eventually get^ the 
value of A^^/VdetA~^. The combination we have access to is totally inde- 
pendent of the amplification factor. As a consequence, the quantity which 
is measurable is the reduced shear field, 

g = 7/(l-«). (75) 

This quantity identifies with 7 only in the limit of very weak lensing (i.e. 
when K <C 1). 



6.4. THE CONSTRUCTION OF THE PROJECTED MASS DENSITY 

The elaboration of methods for reconstructing mass maps from distortion 
fields is not a trivial issue. In a pioneering paper. Kaiser &; Squires (1993) 
showed that this is indeed possible, at least in the weak lensing regime. It 
is indeed not too difficult to show that (in the single lens approximation). 



Vk = - 



( 


82 \ 




V -d2 


di J 


' ( 72 ) 



(76) 



''See Mellier, these proceedings, for a more detailed presentation of the data analysis 
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when K <C 1 and 7^ <C 1. By simple Fourier transforms it is then possible 
to recover k from a distortion map. 

Such a method was further extended in many ways. Bartelmann, Schnei- 
der and his collaborators (Bartelmann et al 1996, Schneider 1995, Seitz &; 
Schneider 1995, Seitz & Schneider 1996, Seitz et al. 1998) have worked 
in detail on the edge effects, the possibility of having adaptive smoothing 
procedures, the use of maximum entropy method etc... This is particularly 
important when structures of different sizes and contrasts are present at 
the same time. Finally Kaiser (Kaiser 1995) exhibited the relation between 
the local convergence and the distortion field g which is valid in all regimes. 
This relation reads 



This is a non-linear and non-local relationship. 

The first reconstruction of a mass map of a galaxy cluster has been 
done on MS1224 by Fahlman et al. (1994). Many other reconstructions 
have now been done or are under preparation (see Mellier 1998 and these 
proceedings) . 

7. The weak lensing as a probe of the Large-Scale Structures 
7.1. THE LARGE-SCALE STRUCTURES 

The idea of probing the large-scale structures with gravitational effects 

is very attractive. The gravitational survey offers indeed a unique way of 
probing the mass concentrations in the Universe since, contrary to galaxy 
survey, it can provide us with mass maps of the Universe that are free 
of any bias. Its interpretations in terms of cosmological parameters would 
then be straightforward and independent on hypothesis on galaxy or cluster 
formation schemes. 

The physical mechanisms are the same in the context of large-scale 
structures and the source for the gravitational effects is the gravitational 
potential given by Eq. (38). It is important to remember that the source 
term of this equation is p(t) (5mass(i) x). The density contrast 5mass(x) is 
usually written in terms of its Fourier transforms. 



The density field is then entirely defined by the statistical properties^ of 
the random variables (5mass(k). At large enough scale the field is (almost) 

^For a detailed introduction to largo-scale structure formation theory and phenomenol- 
ogy see lecture notes of Bertschinger, 1996, and these proceedings. 




(77) 




(78) 
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Gaussian (at least for Gaussian initial conditions which is the case in in- 
flationary scenarios). The amplitude of the fluctuations grows with time in 
linear theory in a known way D^{t). This function is simply proportional 
to the expansion factor for an Einstein-de Sitter Universe. 

The variables are then entirely determined by the power spectrum P{k), 

(<5ma^s(k)5mass(k')) = <5Dirac(k + k') P{k). (79) 

The cosmological model is therefore completely determined by the power 
spectrum, O and A as long as the the dark matter distribution is concerned. 

7.2. THE RELATION BETWEEN THE LOCAL CONVERGENCE AND THE 
LOCAL DENSITY CONTRAST 

The relation between the convergence and the local density contrasts in the 
local universe can be derived easily from Eqs. (38, 52), 

«(7) = ^^^0 / d^. nizs) J dx <^ma.s(x, 7) {l + z). (80) 

In this relation the redshift distribution of the sources in normalized so 
that, 

jdzsn{zs) = l. (81) 

All the distances are expressed in units of c/Hq. The relation (80) is then 
totally dimensionless. 

7.3. THE EFFICIENCY FUNCTION 

It is convenient to deflne the efHciency function, w{z), with 

w{z) = J dzs n{zs) (1 + z) (82) 

On Fig. 7 one can see the shape of the efficiency function for different 
hypothesis for the source distribution. Obviously the further the sources 
are the more numerous the lenses that can be detected are, and the larger 
the effect is. 

7.4. THE AMPLITUDE OF THE CONVERGENCE FLUCTUATIONS 

Prom this equation it is obvious that the amplitude of k is directly propor- 
tional of the density fluctuation amplitude and that the two point corre- 
lation function of the k field is related to the shape of the density power 
spectrum. In the following its amplitude is parameterized with as which 
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Figure 6. Example of reconstructions of projected mass maps. The top panels show the 
initial noise-free n map for either Q. = 1 (left panel) or SI = 0.3 (right panel) with the 
same underlying linear random field and the same rms distortion. The bottom panels 
show the reconstructed k maps with noise included in the shear maps. The maps cover 
a total area of 25 degrees'^. Each pixel has an angular size of 2.5 arcmin^ and averages 
the shear signal expected from deep CCD exposures (about 30 galaxy/arcmin^). The 
sources are assumed to be all at redshift unity and to have a realistic intrinsic ellipticity 
distribution. Such a survey is easily accessible to MEGACAM at CFHT. The precision 
with which the images cam be reconstructed and the striking differences between the two 
cosmological models demonstrate the great interest such a survey would have. 



is the r.m.s. of the density contrast in a sphere of radius 8/i^ Mpc. The 
relation (80) also shows that n depends on the cosmological parameters. 
There is a significant dependence in the expression of the distances but the 
dominant contribution comes from the overall Q,q factor. 

The amplitude of the fluctuations of k depends on the angular scale 
at which the convergence map is filtered. We can introduce the filtered 
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0.5 1 1.5 2 2.5 

redshift 

Figure 7. shape of the efficiency function, w{z) (thick lines), for two different hypothesis 
on the shape of the redshift distribution of the sources (thin fines) . 

convergence kq, with 

i^eil) = J k{j + Y) Wg{j'). (83) 

It is convenient to introduce the Fourier transform of the window function 
W{k). This function is, 

W{k) = 2^, (84) 

where Ji is the Bessel function, in case of a angular top-hat filter. Then 
the filtered convergence reads, 

^eil) = J dxwiz) J ^ ^ ^(k) D^iz) X 

exp [ikrx(z) + ikx.7p(zO)] W [k^9V{z)], (85) 

where the wave vector k has been decomposed in two parts kr and k^ that 
are respectively along the line of sight and perpendicular. The computation 
of the r.m.s. of Kg is analytic in the small angle approximation only. In such 
an approximation we have, 

0<1 kr^k± ^ P{k) Dl{z) P{k±). (86) 

Eventually the variance reads, 

(4) = I dxw'ix) J ^P{k)W\k^eV). (87) 
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For realistic models of the power spectrum (e.g. Baugh and Gaztanaga, 
1995), the numerical result is (Bernardeau et al. 1997), 

/ a \ -(n+2)/2 

{.If' - 0.01 as or (3^) . (88) 

To be noticed is the dependence on the redshift of the sources. This was 
noticed by Villumsen (1996) who pointed out that the VLq dependence is 
roughly given by the value at the redshift of the sources. These results are 
slightly affected by the introduction of the non-linear effects in the shape 
of the power spectrum (Miralda-Escude 1991, Jain & Seljak 1997). 

7.5. THE EXPECTED SIGNAL TO NOISE RATIO 

Arc the effects from large-scale structures measurable? It depends on the 
number density of background objects for which the shape matrices are 
measurable. In current deep galaxy survey the typical mean number density 
of objects is about 50 arcmin"^. The precision of the measured distortion 
at the degree scale is then about, 

for an intrinsic ellipticity of sources of about 0.3. This number is to be com- 
pared with the expected amplitude of the signal coming from the large-scale 
structures, about 1% according to Eq. (88) (see also earlier computations by 
Blandford et al. 1991, Miralda-Escude 1991, Kaiser 1992). This makes such 
detection a priori possible with a signal to noise ratio around 10 (provided 
the instrumental noise can be controlled down to such a low level). 



7.6. SEPARATE MEASUREMENTS OF VL AND (jg 

In Eq. (88) one can see that the amplitude of the fluctuations depend both 
on (Js and on Oq. A question that then arises is whether it is possible 
to separate the amplitude of the power spectrum from the cosmological 
parameters. A simple examination of the equation (80) shows that it should 
be the case, because, for a given value of a^, the density field is more 
strongly evolved into the non-linear regime when is low. 

The consequences of this are two fold. The nonlinearities change the 
angular scale at which the non-linear dynamics starts to amplify the growth 
of structures. This effect was more particularly investigated by Jain &; Seljak 
(1997) who showed that the emergence of the nonlinear regime is apparent 
in the shape of the angular two-point function. This effect is however quite 
subtle since it might reveal difficult to separate from peculiarities in the 
shape of the initial power spectrum. 
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Figure 8. Perturbation Theory results (Bernardeau et al. 1997) obtained for the width 
(left panels) and skewness (right panels) of the probability distribution function of the 
local convergence. The results are plotted as a function of the angular scale in the top 
panels, of fio in the middle ones. The solid lines correspond to source redshifts of 1, and 
the dashed lines to redshifts of 2. In the bottom panels the iso- values of the width and 
the skewness are plotted in the f2o-A plane (for source redshifts of 1). It shows that the 
dependence on fio is slightly degenerate with A. 



The other aspect is that nonhnear effects induce non-Gaussian features 
due to mode couphngs. These effects have been studied extensively in Per- 
turbation Theory. Technically one can write the local density contrast as 
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an expansion with respect to the initial density fluctuations, 

(^ma«s(x) = 4'L(X) + 41s(x) + ■ ■ ■ (90) 

where (5miss(x) is proportional to the initial density field (this is the term 
we have considered so far), (5mass(x) is quadratic, etc. Second order pertur- 
bation theory provides us with the expression of (5mass(x) (there arc many 
references for the perturbation theory calculations, Peebles 1980, Fry 1984, 
Goroff et al. 1986, Bouchet et al. 1993 for the Q, dependence of this result), 

(^ (^ '^lin.(kl)<^lin.(k2) exp[i(ki + k2) • x] X 

"5 ^ ki-k2 ^ 2 (ki-k2)^ 

where 5iin.(k) are the Fourier components of the linear density field. It 
behaves essentially as the square of the linear term, with a non-trivial ge- 
ometric function that contains the non-local effects of gravity. 

Equivalcntly it is possible to expand the local convergence in terms of 
the initial density field, 

«(^) = kW(^)+«;(2)(^) + ... (92) 

The apparition of a non-zero k^^^ induces non-Gaussian effects that can be 
revealed for instance by the computation of the skewness, third moment, 
of kq (Bernardeau et al. 1997), 

{4)={{4^y)+^{{4^y4^)+--- (93) 

The actual dominant term of this expansion is 3 ( '^''^^ ) since the first 

term vanishes for Gaussian initial conditions. For the computation of such 

(2) 

term one should plug in Eq. (80) the expression of SmLs in Eq. (91) and 
do the computations in the small angle approximation (and using specific 
properties of the angular top-hat window function, Bernardeau 1995). 

Eventually perturbation theory gives the following result for a realistic 
power spectrum (Bernardeau et al. 1997), 

^3(^) = 7# = 40 n^'-' z-'-'\ (94) 

The origin of this skewness is relatively simple to understand: as the den- 
sity field enters the non-linear regime the large mass concentrations tend 



(91) 
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Figure 9. Histograms of the values of S3, top-hat fihcr, for = 1 and f2 = 0.3 for a 
5x5 degree survey (thick lines) and a 10 x 10 degree survey (thin lines). 

to acquire a large density contrast in a small volume. This induces rare 
occurrences of large negative convergences. The under-dense regions tend 
on the other hand to occupy a large fraction of the volume, but can induce 
only moderate positive convergences. This mechanism is clearly visible on 
the maps of figure (6). When the mean source redshift grows the skewness 
diminishes since the addition of independent layers of large-scale structures 
tend to dilute the non-Gaussianity. 

What the Eq. (94) demonstrates is that distortion maps can be used to 
determine the cosmic density parameter, provided the redshift distri- 
bution of the sources is well known. The hierarchy exhibited in this relation 
is also a direct consequence of the hypothesis of Gaussian initial conditions. 
Such a hierarchy has been observed for instance in galaxy catalogues (see 
Bouchet et al. 1993 for results in the IRAS galaxy survey). It can be very 
effective in excluding models with non-Gaussian initial conditions (see the 
attempt of Gaztahaga &; Mahonen 1996). To be more precise I present the 
actual histograms of the measured skewness in numerical simulations (Fig. 
9) which clearly demonstrate that the two cosmologies are easily separated. 
One can see that the scatter in S3 is roughly the same in the two cases and 
that the difference in the relative precision is due to the differences in the 
expectation values. 

The validity of Eq. (94) has been confirmed numerically by Gaztanaga 
& Bernardeau (1998), who showed it is valid for scales above a few tens 
of arcmins. A non-zero skewness has also been observed in the numeri- 
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cal experiment of Jain et al. (private communication, in preparation) and 
van Waerbeke et al. (1998). Large angular convergence maps can therefore 
provide new means for constraining fundamental cosmological parameters. 
Numerical results show that in maps of 5 x 5 square degrees it is reasonable 
to expect a precision of a few percent on the normalization and about 5% 
to 10% on the cosmological density parameter depending on the underlying 
cosmological scenario (see Fig. 10). 

7.7. PROSPECTS 

From an observational point of view, the investigation of the large-scale 
structures of the Universe with gravitational lenses is in a very preliminary 
stage. After an early claim by Villumscn (1995), a direct evidence of the 
detection distortion signal of gravitational origin has been reported recently 
by Schneider et al (1997). 

There are at present many studies, either theoretical or numerical, that 
aim to examine all possible systematic errors (Bonnet &: Mellier 1995, 
Kaiser et al. 1995), to optimize the data analysis concepts (such as the pixel 
autocorrelation function by van Waerbeke et al. 1997) and the scientific in- 
terpretations of the resulting mass maps (Bernardeau 1998, Bernardeau 
et al. 1997, Seljak 1997, van Waerbeke et al. 1998). A few observational 
surveys are now emerging, 

— the ESO key program jointly done by MPA and lAP; 

— the DESCART project, part of the scientific program of the wide field 
CCD camera to be installed at the CFHT. 
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